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The Sunyaev-Zel'dovich effect due to hyper-starburst galaxy winds 
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ABSTRACT 

We construct a simple, spherical blastwave model to estimate the pressure structure of the 
intergalactic medium surrounding hyper-starburst galaxies, and argue that the effects of inter- 
action with star- forming galaxy winds may be ap proximated at early t imes by an adiabatically 
expand ing, self-similar 'bubble' as described by IWeaver et alJ (Il977h and lOstriker & McKed 
(Il988l) . This model is used to make observational predictions for the thermal Sunyaev- 
Zel'dovich effect in the shocked bubble plasma. Radiative cooling losses are explored, and 
it is found that bremsstrahlung will limit the epoch of adiabatic expansion to 10^-10'* years: 
comparable to total hyper-starburst lifetimes. Prospects for making a first Sunyaev-Zel'dovich 
detection of galaxy wind bubbles using the Atacama Large Millimeter Array are examined for 
a number of active hyper-starburst sources in the literature. 

Key words: shock waves — galaxies: high-redshift — galaxies: starburst — intergalactic 
medium — cosmic microwave background 



1 INTRODUCTION 

After the epoch of reionization, free electrons in the intergalac- 
tic medium (IGM) interact with photons in the cosmic microwave 
background radiation (CMB) via Thomp son scattering. Known 
as the Suny aev-Zel' dovich (SZ) effect (e.g . ISunvaev & Zel'dovichl 
I1970II1972I : see also lCarlstrom et al.ll200a for a review), the result 
is a characteristic spectral distortion of the CMB. The effect can be 
split into two contributions: that due to the thermal motion of free 
electrons which causes an overall increase in the apparent tempera- 
ture of CMB photons (called the thermal SZ or tSZ effect), and tha t 
due bulk kinetic motion (the kSZ effect; see lCarlstrom et al.l2002h . 
The tSZ effect, in particular, shows promise as a means of 
observing the hot plasma surrounding galaxy clusters and other 
peaks in the matter power spectrum, and recent ye ars have seen a 



Solomon & Vanden Bout 2005'; IWang et al.' 2008; ICasey et al.l 
20091; iMartinez-Sansigre et al.l |2009| ; iRiechers eTal] |2009|; 



Wagg et alj 120091) . Other recent results suggest that the galaxy 



steady growth in both theoretical predictions (e.g. Moodley et al 
_2009; Scannapieco, Thacker & Couchman 2008; Roncarelli et al 



wind outfl ows frequently ob served in star-forming galaxy spectra 
(see, e.g., IVeilleux. Cecil & Bland-Hawthorn 2005 for a review) 
are ub iquitous wherever there is star format i on activity to drive 
them JGrimes et alj|2009l ; iRubin et al.ll201(]| ; IWeiner et al.ll2009l; 
iKrug et al.ll201(]h . with wind masses and velocities increasing with 
M, . The winds associated with hyper-starbursts at high redshift 
will be extremely energetic, leading to violent interaction with the 
plasma of the surrounding intergalactic medium (IGM). 

The following question naturally arises: can we de- 
tect the tSZ due to the wind-IGM interaction for these 
events? A s imilar question has been posed for active galax- 
ies ( e.g. J Yamada. Sugivama & Silk Il999l; IPlatania et al.l 



20070 and actual measurements for massive galaxy clus ters (e.g. 
Halverson et"Zll2009l ; Isasu et al.ll2oIol ; lNord et al.ll2009l) . It is in- 



2OO2I; IChatteriee & Kosowskvl |2007|: 



Scannapieco et al. 2001 



teresting to consider what other astrophysical objec ts may be ob- 
serva ble in the near future with the tSZ effect (e.g. lYamada et al.l 
I20T0I) . 

In recent years an increasing number of high redshift 'hyper- 
starbursts' (which we define as objects with an estimated star 
formation rate of M, > 10^ Moyr"'^) have been reported in the 
literature, often but not always identified with submillimetre galax- 
ies (SMGs), and are believed to represent a s ignificant fraction 
of the total star formation at these epochs (e.g. iBlain et al.ll2002l ; 
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Chatteriee et al.ll2008l ; lYamada et al.ll2010|) and for th e combined 
signa l due to star-forming winds (e.g. iMaiumdar. N ath & Chiba| 
l200ll ; lBabich & Loebll2007h . but individual hyper-starburst objects 
now also raise interesting possibilities. The temperature increase 
of C MB photons due to th e tSZ effect is redshift independent 
(e.g. ICarlstrom et al.l 120021) . which potentially opens a new 
observational window on extreme star formation at high redshift. 

We examine this question by constructing a simplified adia- 
batic model of the wind-IGM interaction (Section[2]l, and then cal- 
culate the tSZ effect according to such a model for a variety of pos- 
sible hyper-starburst targets in the literature (Sections[3]&|4ll along 
with observational prospect s for the upcoming Atacama Large M il- 
limeter Array (ALMA: see lBrown. Wild & Cunninghanill2004t) fa- 
cility. In Sections[5]we test the adiabatic expansion assumption by 
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calculating radiative cooling timescales for the hot plasma, and end 
in Section|6]with a discussion and conclusions. 



2 GALAXY WIND INJECTION BUBBLE MODEL 

As demonstrated bv lSamui. Subramanian & Srianandl OOOSO . who 
present a detailed semi-analytic analysis of galactic outflows, the 
galaxy wind-IGM interaction can be approximately described at 
early times as a spherically expanding, self-similar blastwave or 
'bubble'. T he model used to des cribe this bubble is based on that 
presented in lWeaver et al.l ( Il977l) . which describes the adiabatic in- 
teraction of strong stellar winds with the interstellar medium. This 
is related to the original Sedov-Taylor blastwave sol ution (see, e.g., 
iLandau & Lifshitzlll959l ; lOstriker & McKeel Il988l) but describes 
the shock evolution under continuous rather than instantaneous en- 
ergy injection, making it being better suited to cases with ongo- 
ing starburst activity. The I Weaver et al.l ( 119771) model retains self- 
similarity by not introducing a characteristic timescale to the stel- 
lar wind bubble, and it is this property that allows it to be applied 
equally well to the galaxy wind-IGM interaction despite the greater 
energies of this regime. 

At time f = star formation processes initiate the expul- 
sion of a galaxy wind, which we approximate as having constant 
mass outflow rate Afw leading to a mechanical energy injection 
rate given by 



Lw = M„vl/2 



(1) 



where v^ is the net wind outflow velocity. All quantities 
will be expressed in relation to a fiducial value Ww ~ 10'^ 
kms~^, motivated by observations of intermediate-to- high redshift 
wind s for objects with extreme star formation rates jRubin et al 



lOldlWeiner et al.l2009l:ICapak et alj2008l ; lBland-HawthoiTi et al 



2OO7I : 'Adelberger et alj l2003h . Following lAguirre et all 2001 ) 
Sprin gel & HernquistI l l2003h . and lDalla Vecchia & Schavd hoosi) 
we relate the wind mass outflow M^ to the star formation rate il/* 
using the dimensionless wind mass loading parameter r; (adopting 
a fiducial value of 1 in all calculations), so that Afw = ryA/,. The 
total energy injected by the wind into the system is thus 



E = L^t = rjA'Lvlt/2. 



(2) 



The mechanical luminosity Lw can be related to the total energy 
available from core collapse supemovae via the wind energy frac- 
tion 



/w = 77«w/(2esN), 



(3) 



whe re 6sn is the kinetic e nergy injected per Mp of stars f ormed 
(e.g. lDalla Vecchia & Scha ve 2008). Assuming a lChabriej | |2003|) 
initial mass function, these authors adopt a value esN — 1.8 x 
lO^^ergMg'^ : for our fiducial values of v,, and 77 this implies a wind 
energy fraction of /„ ~ 0.55. This is slight l y larg er than the value 
of 0.4 adopted bv lOalla Vecchia & Schavd ( 120081) for more mod- 
erately star-forming systems. It should be cautioned that the ex- 
act value of /w appropriate for supernovae in high-redshift hyper- 
starburst regions is unclear, and so all results will be presented with 
scaling relations for the parameters v^ and 77. 

The interaction of the wind with the IGM can be split into four 
physically distinct zones at increasing radial distance r from the 
bubble centre: a) the free-streaming wind immediately on exit from 
the hyper-starburst region; b) a region of shocked galaxy wind; c) 
a shell of shocked IGM plasma; d) the ambient IGM surrounding 
the hyper-starburst host galaxy l lWeaver et alj 119771 ; ISamui et al.l 



120081) . We model the unshocked IGM in zone d) as ionized hy- 
drogen of homogeneous density pd- In order to estimate a value 
for this density in a hyper-starburst galaxy, we make use of mea- 
surements for t he clustering bias feq of high-redshift SDSS-DR5 
quasars given in lShen et alj ( 120091) : large molecular gas reservoirs 
and hyper-starburst events are often associated with quasi-stellar 
objects (QSOs ), radio galaxies, or other indicat ions of active galac- 
tic nuclei (e.g. lSolomon & Vanden Boutll2005l ; see also Section |4] 
Table [1). These authors find a bias of 6q = 12.96 ± 2.09 for the 
sample of 1788 objects in their highest redshift bin: 3.5 < 2 < 5.0 
(median z = 3.84), and we thus choose a fiducial modelling value 
of 6q — 13. We also assume that the IGM in the region immedi- 
ately surrounding the starburst is biased in the same way, giving 



Pd ~ Pcritrib(l + Z) (1 + BqS), 



(4) 



where pcrit and fib are the current-epoch critical density and 
baryon density parameter respectively, and 5 is the dimensionless 
density perturbation. We approximate this perturbation in the star 
formation environment as <5 ~ 180, the mean overall matter density 
perturbation for a collapsed halo. 

The outer radius of an adiabatically expanding (non-radiative), 
spherical blastwave under constant energy injection is given by 



i?2(t)=/?(LwtVpd)' 



(5) 



whe re P = 0.8828 is a const ant estimated via numerical calcula- 
tion JOstriker & McKeel 1988h . Equation ([5} is the only dimension- 
ally correct combination of the system variables in this time and 
distance scale-free problem, and its power law form is therefore 
required; this self-similar evolution will cease only after radiative 
losses become significant (see Secti on [5). Using a flat. W MAP 5- 
year best-fitting ACDM cos mology (IKomatsu et al.ll2009l) and the 
value fib = 0.044 ± 0.01 dSanchez et al.ll2009h . we may rewrite 
equation (O in terms of fiducial values for a hyper-starburst shock- 
wave system as follows: 



i?2(i) 



/3 



■qM^vlt^ 



= 28.86 r?s 
t 



2nl,p,,it{l + z)^l + bQS) 



(6) 



W Vl03kms-V 



lO^yr 



lO^Moyr-V VlOakms^^ 

_ 3 

(1 + 2:) 5 kpc. 



(7) 



It can be seen that bubble radius will quickly extend significantly 
beyond the star forming regions into the surrounding IGM. 

The pressure and den sity structure within the bubble is als o 
approximately self-similar dWeaver et al.ll977l ; ISamui et alj2008h . 
and can therefore be described in terms of the dimensionless radial 
parameter ^ — r/R2{t). At the very edge of the shock (^ — 1), the 
gas quantities are determined by the Rankine-Hugoniot conditions 
for a strong shock boundary: 



P2 = pd{j + 1)/(7 - 1), P2 = 2pdm/{j + 1), 



(8) 



where R2 = d_R2/di, and for the monatomic gas w e have adi- 
abatic index 7 — 5/3 (see lUandau & Lifshitjl 19591) . Using this 
value gives an edge-of-shock density of p2 = 4: pd and pres- 
sure P 2 = 3pdi?2/4. Within the outer shock radius. IWeaver et alj 
( 119771) find that the density p(^) decreases inwards and drops sud- 
denly to zero at the contact surface Re = 0.86i?2, the interface 
between the shocked IGM in zone c) and the shocked wind in zone 
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Figure 1. Pressure (solid line) and Compton y parameter (dashed line) in 
the model as a function of distance from the bubble centre, with labels for 
the four distinct physical zones a)-d). 



b). The pressure at this surface is P{^ = 0.86) — 0.59/9di?2, and 
these scalings hold while the expansion remains adiabatic. 

For a perfect gas the tSZ effect may be estimated by inte- 
grating the pressure along a line-of-sight path through the bub- 
ble; a model for the pressure structure throughout must there- 
fore be adopted. An exact model ma y be calculated numerically 
jWeaver et al.lll977l : lDokuchae\il2002n but ultimately the accuracy 
of such a model will rely on the physical conditions matching those 
described above (e.g. uniform ambient IGM density pd or constant 
input star-formation power L„), which is unlikely. Instead, for the 
purposes of estimating an approximate observability of the tSZ ef- 
fect we adopt a simple model that displays the correct behaviour at 
the boundaries and makes simplifying interpolations between them: 



PiO 



Pb 

Pb + 



< C < 0.86 
^^^^ "> (C - 0.86) 0.86 < C < 1 ^^'' 



where Ph = 0.59pd^2- The model, plotted in Figure[T] is isobaric 
for r < Rc\ this b ehaviour is argued from simple principles by 
IWeaver et alj ( 11977!) and can also be see n in the complex analytic 
solutions provided bv lPokuchaevI d2002r) . This pressure model can 
then be used to estimate the tSZ effect for a hyper- starburst wind 
bubble. 



3 THERMAL SZ EFFECT IN THE BUBBLE 

The tSZ effect describes the spectral variation of the specific in- 
tensity of the CMB photons caused by their pa ssing through a 
non-r elativistic. Maxwell-distributed electron gas l lCarlstrom et al.l 
l2002h . It takes the form 



AT, 



a^expa; 



y 



(10) 



exps — 1 \^tanha;/2 

with X = hv/kBTcMB (using Tcmb = 2.725 K at 2 = 0), and 
where the Compton y parameter is defined via the path integral 

keT^il) 



y = 



I 



-(TTne(Z)dZ 



(11) 



along the observer line-of-sight. We do not include a relativis- 
tic correction to the tSZ effect (e.g. IChallinor & Lasenbvlll998l ; 
lltoh, Kohvama & Nozawalll998r) for the hot Shockwave electrons 
(see Section|4l(. 



For a perfect gas such as the extremely rarefied IGM, equation 
dllb becomes 



ksPil) 



ajdl. 



(12) 



Using the fiducial values described, and the pressure model de- 
scribed by equation ([9]l, this gives a peak value of 



ypeA = 2.456 X lO"" 775 



M, 



lO'^yr 



10^ Moyr" 
(1 + ^)* 



-Y (^ 1 

W Vl03kms-V 



(13) 



along the line-of-sight through the centre of the bubble. Note that, 
due to the dependence of P upon pd, the strength of the tSZ in this 
model increases with redshift (although the bremsstrahlung cooling 
rate will also: see Section[5]l. The y at any distance from the bubble 
centre can be calculated using the expression for pressure given in 
equation ([9]l, and is shown in Figure[T] 



4 POTENTIAL HYPER-STARBURST TARGETS 

In Table[T]we present a limited selection of hyper- starburst objects 
from the literature that may make suitable targets for tSZ observa- 
tion, along with predictions for observable properties based on the 
simple bubble model. In all these calculations we use the published 
estimate of the SFR, take iiw = 10"^ kms~^ and 77 = 1 (correspond- 
ing to /w — 0.55), and employ a flat, WMAP 5-year cosmology as 
in Section|2] Where an estimate of the age t of the hyper-starburst 
exists (e.g. from SED fits) that value is chosen; otherwise, the de- 
pletion timescale for the CO-detected molecular gas reservoir is 
used as a proxy (since on average a hyper-starburst will be observed 
midway through this depletion). The observational angular radius 
62 corresponding to R2 is given for each bubble model, and the 
central line-of-sight Compton j/peak is tabulated. In all cases this is 
seen to be significantly greater than the cosmolo gical expectation 
value of (y) = (1.19±0.32) x 10"^ calculated bv lRoncareUi etall 
l l2007h . 

Using the ALMA sensitivity calculator (made available on 
the European Southern Observatory web pageo we also estimate 
SA^IJ , defined as the surface brightness signal-to-noise that may 
be obtained at j/pcak using ALMA within an exposure time of 24 
hours. We calculate this value assuming continuum observation 
centred on ;/ = 400 GHz (within ALMA Band 8), and an antenna 
configuration giving a beam with FWHM = 0.75 arcsec. In calculat- 
ing SN24 each bubble is first convolved with this simulated Gaus- 
sian beam to estimate the observable brightness temperature excess 
at this resolution. The sensitivity at the continuum band centred on 
f is calculated using the online observing tool, taking the declina- 
tion (and hence typical atmospheric column density) of each object 
into account. For those targets with a sky position that cannot be 
observed using the Southern Hemisphere ALMA facility (5 ~ 30° 
and above) we instead use an equatorial location for comparison. 

To visualize these results we simulate an ALMA image of 
one of the bub bles visible from the Southern Hemisphere, SMM 
J222 1 74-1-00 15 JNeri et alj|2003l ; ISolomon & Vanden Boulll2005l) . 
at the same v and beamsize but for an extended exposure time of 

5 X 24 hours: this bubble can be seen in Figure |2] The map also 
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Table 1. Model predictions for a selection of reported hyper-starbursts, based on published redshifts and star formation rate estimates (see Sections|4]&|D- 



Galaxy target ' 



Redshift Type SFR t R2 82 j/pcak 

z MQyr^-"^ lO^yr kpc arcsec 10^^ 



RG J123649.66+620738.0'' 

SMMJ222174+0015'' 

4C41.17R/B^ 

SMMJ154137+6630.59 

COSMOS J100054+023436'"'' 

BR 1202-0725-' 

LESS J033229.4-275619'= 

SDSS J092721.82+200123.7' 

SDSS J114826.64+525 150.3" 



2.32 
3.10 
3.80 
3.93t 
4.55 
4.69 
4.76 
5.77 
6.42 



Radio 
SMG 
Radio 
SMG 
SMG 
QSO 
SMG 
QSO 
QSO 



3800= 
1800= 
3000= 

> 10 000" 
2900*" 
9000= 

~ 1000*= 
3200™ 
3300° 



IT 
1.7= 
1.1 = 
It 
0.7'^ 
0.9= 

0.5"" 
IP 



lO^yr lO^'yr 



18.4 
19.1 
14.9 
17.6 
10.3 
14.8 
23.2 
7.63 
11.1 



2.21 
2.47 
2.06 
2.46 
1.54 
2.27 
3.53 
1.28 
1.97 



2.32 
1.71 
3.06 
6.63 
3.91 
7.55 
1.52 
5.63 
5.67 



1.19 
0.88 
1.56 
3.40 
1.95 
3.97 
0.82 
2.42 
2.89 



971 
410 
219 
200 
124 
112 
106 
56.3 
38.8 



36.0 
11.6 
8.73 
10.5 
6.14 
6.37 
2.16 
3.50 
1.93 



on or spec tra l confirmat i on); iChapman et al.l J2004h: ICasev et alj 2005); JNerie^al^ 

: et al' l'2005"); s'Wagg et al.' (2009); ''Capak et al. l'2008"); ''Schinnerer et al." 200f ); - lOmont et aj 



(References for CO/FIR-emis si on detection 
ISolomon & Van den Bout ('2005^; ^De Breuck f 

ICoDDinetalf l2009); 'CariUi et al. (2007); '"Wanget al. (2008); "Bertoldi et al. (2003); "^Riechers et_aL, 
tphotometric redshift estimate. -^Unknown, fiducial value used (see equations ItI& I 1 3 I for scalings). 



20 08|); 40mont et alj 
J^ .Walter et al. c2Q0l . 



2003) 



19961) 



includes contamination due to warm dust emission in the hyper- 
starburst SMG, for which we assume dust grains at 40 K with emis - 
sivity £1, oc u^ and spectral index 13 = 1.5 (e.g. lBlain et alj|2002h . 
We model the emitting surface brightness distribution as a uniform 
disc of diameter 2 kpc in the centre of the bubble, comparable 
with r ecent observations of compact hyper-star bursts dCasev et al.l 
I2OO9I ; IWalter et al.ll2009l ; iMaiolino et al.ll2007h . Taking the value 
of 6 3 ± 1.3 mJy for the SMM J222 174+00 15 flux density at 850 
^m d Solomon & Vanden Boutll 20051) we estimate a flux density of 
9.7 mJy at 400 GHz, and plot contours of the convolved bright- 
ness temperature distribution in the central region of Figure |2] The 
brightness temperature greyscale is here removed subject to a sim- 
ple cut of ATt < 0.4 mK showing that, thanks to the angular res- 
olution of ALMA, the extended bubble is discernible beyond the 
strong-but-localized dust emission at its centre. In practice, a more 
sophisticated removal of SMG dust contamination will be aided by 
the possibility of ultra-high resolution mapping of the bright cen- 
tral dust (ALMA can reach FWHM = 0.011 arcsec at 400 GHz), 
and by the characteristic signature of the tSZ effect on CMB pho- 
tons: a maximum spectral exces s/decrement arou nd 385/144 GHz 
and null crossover at 218 GHz dCarlstrom et alj2002i) . Multiple ob- 
servations will require further telescope time, however. 

In Section [3] it was noted that relativistic corrections to the 
tSZ are not included in the determination of y. Using equation 
l[8j, the temperature of the gas in the shock boundary is T2 = 
3(me + mp)7?2/(4fcB), which implies atemperature of r2 > lO^K 
for all the objects listed in Table [T] Calculating the reduction in 
tSZ intensity using the first and second order relativistic terms de- 
scribed by Challinor & Lasenby ( 199^ yields corrections of ~ 3% 
for the hyper-starbursts listed. The effect is therefore negligible for 
the current model, but relevant for more detailed calculations. 




X [arcsec] 

Figure 2. Simulated ALMA excess brightness temperature map at 400 
GHz, relative to the CMB, for SMM J222174+0015 (as described in Sec- 
tion|4j. The simulated exposure time is 5 X 24 hours, and contours in the 
greyscale-shaded region are placed at Icr intervals around ATb = 0. We 
remove the greyscale map in the central, SMG dust-contaminated region 
subject to a threshold of AT;, < 0.4 mK, but contours are here plotted at 
intervals of 100 mK to illustrate this bright central source. The hatched area 
in the upper left illustrates the simulated beam of FWHM =0.75 arcsec. 



5 RADIATIVE COOLING IN THE INJECTION BUBBLE 

The discussion so far has assumed adiabatic bubble expansion, but 
depends upon losses due to radiative cooling being small within 
the hyper-starburst age t. At the shock boundary temperatures of 
T2 > lO^K the gas will be fully ionized but not sufficiently en- 
ergetic for electron-positron annihilation to contribute to radiative 
processes: thermal bremsstrahlung and inverse-Compton cooling 
by CMB photons will dominate radiative losses (this latter being 
the net cooling effect of the tSZ itself). The radiative cooling will 
be greatest at the shock boundary R2, where the combination of 



plasma temperature T2 and density p2 are greatest, and so we cal- 
culate timescales for radiative losses in this region to explore the 
validity of the adiabatic expansion approximation. 

We consider first the effect of inverse-Compton scattering of 
CMB photons. In the relevant non-relativistic limit, the net energy 
loss rate for a single electron in the shock plasma is given by 



dE 16 . (vl) 



(14) 



(see, e. g. , lLongaiii 1 992h where aj is the Thomson scattering cross- 
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section, (jsb is the Stefan-Boltzmann constant, Tcmb is the CMB 
photon temperature at that epoch, and Ve is the velocity of the elec- 
tron. This gives an exponential decay law with an inverse-Compton 
cooling timescale of 

TIC = (3meC^) / [32aTasBTcMBiz)] (15) 

at the epoch at which the shock is observed. Values of ric for the 
selection of targets discussed in Section|4]can also be seen in Table 
[T] and in all cases can be seen to be large in comparison to the 
typical hyper- starburst lifetime of ~10^ years. 

Cooling due to thermal bremsstrahlung near the shock bound- 
ary i?2 may be more significant. A thermal hydrogen electron in 
the shock boundary plasma will lose energy at a rate 



-^ = 1.435 X 10'*°npgVT2, 



(16) 



where rip is the volume number density of pro tons in the pla sma 
and g ~ 1.2 is taken as a typical Gaunt factor l lLongaig|l992n . In 
Table [T] we use this relation to give timescales tb over which the 
shock boundary energy would be reduced to 37% of its instanta- 
neous adiabatic value by bremsstrahlung alone, as an indicator of 
the importance of radiative cooling. These timescales are generally 
longer than the typical lifetimes of hyper-starbursts, but not always: 
bremsstrahlung will clearly dominate over inverse-Compton losses, 
and at higher redshift (where Hp is greater) will limit the epoch of 
pure adiabatic evolution to within t < 10^ years, comparable with 
the duration of hyper- starburst activity. However, these radiative 
losses are not total: the adiabatic model can be viewed as giving 
an approximate upper bound on the expected range of j/pcak. This 
increases the importance of bubble age in the tSZ signal beyond 
the otherwise weak y oc t~^'^, and might be used to provide lower 
bounds on t should a tSZ bubble not be detected. 



6 DISCUSSION & CONCLUSIONS 

We have presented an adiabatic model of wind-driven bubble 
Shockwaves around hyper-starburst galaxies and made predictions 
for the strength of the tSZ effect due to such objects. Assuming this 
model to be a reasonable approximation we have shown that such 
bubbles are detectable using ALMA, albeit requiring integration 
times of multiple days (in two or, preferably, three bands). We have 
shown that bremsstrahlung will dominate the radiative cooling of 
the shocked gas, and found that this will cause the evolution of the 
shock to become significantly non-adiabatic after a time compara- 
ble to the typical total duration of hyper-starbursts themselves. 

In order to construct the model a number of simplifications 
have been made. Among these, the assumption of spherical sym- 
metry and the homogeneous density pa are least likely to be valid. 
The former assumption may be thrown into question for systems 
with active galactic nuclei, as are often seen for hyper-starbursts, or 
with star formation occurring in an exten ded disc rather than a cen- 
tral core: lDalla Vecchia & Schavd ( 120081) found a bi-polar pressure 
distribution when modelling outflows in massive disc-like forma- 
tion, once winds were fully coupled to the gas dynamics of the sys- 
tem. A bi-polar outflow would alter the pressure distribution of the 
shocked gas and might reduce the strength of the tSZ effect along 
certain lines of sight. 

However, the star typical star formation rates in the mas- 
sive systems simulated by these authors were small when com- 
pared to the figures in Table [T] and it might be expected that 
the greater mechanical energy input from hyper-starbursts will al- 
ter the wind-IGM interaction significantlv. lOalla Vecchia & Schavd 



( 120081) found interesting qualitative differences in the behaviour 
of outflows for their simulated disc systems of total halo mass 
lO^°h^^MQ and lO^^/i^^M© that produced early star formation 
rates of ~ 0.05 M0yr^^ and ~ 10 M0yr^^ respectively. An in- 
crease of a further two orders of magnitude in star formation rate 
may cause sufficient increases in wind ram pressure as to qualita- 
tively change the outflow beha viour once more. In a related study 
IChatteriee & Kosowsk^ ( |2007|) and lChatteriee et alj |2008!) found 
that quasar-driven bubbles in simulated halos provided a reasonable 
match to a simple spherical model, although again at lower energies 
and typically after longer timescales than we consider here. Unfor- 
tunately there is insufficient knowledge and modelling of the wind- 
interstellar medium interaction in these powerful systems to do 
more than speculate. The typical morphologies and environments 
of the star-forming gas regions prior to the onset of strong winds are 
also poorly known for these hyper-starbursts, which further com- 
plicates accurate simulation. More detailed theoretical modelling 
combined with future precision observations of the tSZ might help 
shed light on exactly these questions. 

Another simplification, the assumption of a homogeneous 
IGM, may misrepresent the evolution of the bubble with time: for a 
more realistic density profile the bubble will initially expand more 
slowly through the dense gas at the halo centre, and then more 
rapidly at later times as pd{R2) decreases. This is an area in which 
the modelling could be improved, particularly if aided by more de- 
tailed knowledge of the gas density around hyper-starbursts. A re- 
lated simplification is the model's neglect of the gravitational de- 
celeration of the Shockwave due to the host halo potential. This 
effect will generally represent a small correction to the dynamics 
of adiabatic expansion and radiative losses, which can be seen by 
comparing the galaxy wind velocities to the escape velocities of 
most galactic haloes (entrained clouds, however, might impart a 
velocity kick that prevents their escape). These effects should be 
considered in more accurate modelling, perhaps using a simpli- 
fied model of a galaxy halo potentia l and mass entrainment (e.g. 
iMurrav, Ouataert & Thompsodl2005h . 

Important uncertainties also arise due to the assumed val- 
ues of the starburst age i, and the fiducial parameters ?)„ and 
AI, (or equivalently /„). It should be emphasized that the fidu- 
cial predictions may be adapted as appropriate using the simple 
scaling relations presented, which may be necessary as we gain 
a clearer picture of hyper-starburst physics. As an example, using 
y oc /w from equation J13[l an d adopting the value /„ — 0.4 
of' Dalla Vecchia & Schavd ( 1200 8h reduces the tSZ signals listed in 
Table[T]by ~ 20%. However, a drastic reduction in /„ (appropri- 
ate if the kinetic energy from hyper-starburst supemovae is mostly 
lost radiatively) mi ght yet reduce y t o the c osmological confusion 
limit calculated by iRoncarelli et al.l ( 120071) . Conversely, observa- 
tions of the tSZ around hyper-starbursts could in this way, even if 
not detected, provide interesting constraints on /„ and thus super- 
nova wind feedback in these objects, whi ch may contribute much 
of the total star-formation at high redshift dBlain et al.l2002h . They 
might also constrain the feedback energy injected by the active 
galactic nuclei comm only hosted in hyper-starburst systems (e.g. 
IChatteriee et al.l2008r) . although more careful modelling will be re- 
quired to distinguish the two signals. 

To conclude, wind-driven bubbles around hyper-starburst ob- 
jects may be detectable via the tSZ in the near future using submil- 
limetre instruments such as ALMA or its successors. Despite many 
simplifying assumptions in the model, the enormous mechanical 
energy of the wind-IGM interaction will ensure that a system ap- 
proximating the simple, self-similar blastwave must persist while 
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radiative losses are small. According to such a model, such ob- 
jects as those listed in Table [T] (with estimated j/peak > 2 x 10"'') 
might be securely detected with ALMA within days. The number 
of known hyper-starbursts in the Southern Hemisphere must be ex- 
pected to increase with time, and one of these may present an ex- 
cellent opportunity for a first bubble detection with ALMA. 
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